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V  SUMMARY 


->The  results  of  calculations  of  the  flow  over  slender  circular  cones  at  large  angles 
of  Incidence  are  presented.  The  vortices  which  form  above  the  cone  as  a  result  of 
boundary- layer  separation  are  represented  by  two  Invlscld  models:  as  Isolated  line- 
vortiees  and  as  vortex  sheets  with  llne-vortex  cores.  For  each  model,  an  essentially  new 
family  of  solutions  Is  revealed,  in  which  the  flow  remains  asymmetric  even  when  the 
separation  lines  are  specified  to  lie  sysmetrically  about  the  Incidence  plane .  The  side 
force  which  is  predicted  by  this  second  family  of  solutions  is  very  much  larger  than  the 
side  force  which  can  be  obtained  by  specifying  asymmetric  separation  lines  in  the  first 
family  of  solutions,  and  corresponds  to  the  maximum  side-force  levels  found  experimentally. 
It  is  concluded  that  an  essentially  invlscld  mechanism  for  the  generation  of  side  force 
at  zero  yaw  has  been  discovered. 


LIST  OF  SYMBOLS 


Cp  pressure  coefficient 

Cy  side-force  coefficient 


s 


normal-force  coefficient 


i 


r  radius 


0  freestream  speed 

x,y,z  Cartesian  body  axis  system  (with  perturbation  velocities  u,  v,  w) 

z  complex  coordinate  (»  y  +  i*) 

a  angle  of  incidence 

r  circulation 

6  semi-apex  angle  of  cone 

8  angular  coordinate 


subscripts 

1,2  right  hand,  left  hand 

1,2, 3, 4,5  complex  planes  in  transformation  sequence 

L,R  left,  right  (vortices) 


1  INTRODUCTION 

This  paper  is  concerned  with  the  out-of-plane  force,  or  the  side  force  at  zero  yaw, 
which  arises  on  slender  bodies  at  large  angles  of  incidence.  Space  does  not  permit  a 
o.  detailed  review  of  the  phenomenon,  but  reference  may  be  made  to  the  recent  account  by 

-«  Ericsson  and  Reding1 ,  and  to  the  papers  they  cite.  There  is  at  present  no  general  agree- 

_  ment  about  the  mechanism  responsible  for  the  creation  of  out-of-plane  forces,  and,  indeed, 

more  than  one  mechanism  may  be  involved  in  the  whole  range  of  body  shapes,  angles  of 
u  incidence,  Mach  number  and  Reynolds  number  at  which  this  phenomenon  is  observed.  However, 

4j  there  is  agreement  that  the  phenomenon  is  associated  with  flow  separation  from  these 

bodies,  and  that  the  detailed  shape  of  the  nose  of  the  body  has  a  strong  influence.  In 
P  general  terms,  at  a  sufficiently  large  angle  of  incidence,  the  flow  on  the  body  surface 

diverges  from  an  attachment  line  on  the  windward  side  and  separates  along  two  lines  on  the 
sides  of  the  body.  The  shear  layers  which  spring  from  the  separation  lines  roll  up  into 
a  pair  of  contra-rotating  vortices  which  lie  above  the  body. 
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We  shall  only  concern  ourselves  with  pointed  bodies  of  revolution  at  angles  of 
incidence  for  which  the  separation  lines  extend  forward  to  the  pointed  nose.  In  these 
circumstances,  the  vortices  form  at  the  nose,  and  v..eir  initial  development  is  the  same 
as  it  would  be  on  a  circular  cone.  For  a  certain  range  of  incidence,  the  flow  is  always 
found  to  be  symmetric,  but  at  larger  incidences,  as.mmetric  flows  are  observed.  In 
particular,  the  two  separation  lines  on  the  surface  are  not  symmetrically  placed  about 
the  incidence  plane  and  the  two  vortices  also  lie  in  asymmetric  positions.  The  question 
arises  whether  it  is  the  asymmetry  in  the  separation  line  positions  that  causes  the 
asymmetry  in  the  vortices,  or  the  asymmetry  in  the  vortices  that  causes  the  asymmetry  in 
the  separation  lines.  There  is,  of  course,  a  sense  in  which  it  is  inappropriate  to  speak 
of  cause  and  effect,  but  we  believe  the  evidence  presented  here  will  show  that  there  is 
an  inherent  asymmetry  in  the  inviscid  flow  to  which  the  viscous  flow,  and  hence  the 
separation  positions,  must  adjust.  In  other  words,  the  essential  mechanism  responsible 
for  the  out -of -plane  forces  is  an  inviscid  one. 

Evidence  for  this  will  be  given  in  the  form  of  results  of  calculations  using  two 
simple  inviscid  models  of  the  flow  past  a  circular  cone  at  incidence,  together  with  3ome 
recent  measurements  made  at  RAE  by  Mundell2.  Since  the  theoretical  models  are  inviscid, 
the  positions  of  the  separation  lines  are  specified  at  the  start  of  the  calculation. 

In  principle,  the  separation  positions  could  also  be  predicted  for  laminar  flows,  using 
the  method  described  in  an  earlier  paper  by  Fiddes3  for  symmetric  flows.  However,  we 
have  exploited  the  freedom  to  specify  the  separation  positions. 

In  the  first  place,  we  have  started  with  a  previously  known  solution  for  the 
symmetric  flow  past  a  circular  cone,  in  which  the  separation  position  has  been  specified 
symmetrically.  This  solution  has  then  been  modified  progressively  by  introducing  a 
gradually  increasing  asymmetry  in  the  separation  lines.  This  gi-  .s  rise  to  a  moderate 
degree  of  asymmetry  in  the  vortices,  as  regards  their  positions  and  circulations,  and1*! 
also  produces  a  side  force.  However,  even  for  a  large  asymmetry  in  the  positions  of 
the  separation  lines,  the  side  force  remains  a  small  proportion- of  the  normal  force, 
whereas  measured  values  of  the  side  force  are  known  to  approach  the  magnitude  of  the 
normal  force.  This  leads  us  to  believe  that  the  phenomenon  is  not  being  driven  by  a 
boundary  layer  mechanism,  such  as  differential  transition. 

In  the  second  place,  we  have  specified  symmetric  separation  line  positions  on  the 
circular  cone  and  found  a  family  of  solutions  in  which  the  vortices  are  asymmetric. 

Results  have  already  been  published4  for  the  simpler  of  the  two  mathematical  models 
considered  in  which  the  vortices  are  represented  by  line  vortices  joined  to  the  separa¬ 
tion  lines  by  cuts  (a  model  first  used  for  the  symmetrical  flow  past  bodies  by  Brysor.s)  . 

New  results  have  been  obtained  for  the  more  realistic  model,  described  for  symmetric 
flows  by  Fiddes3,  where  the  shear  layers  are  represented  by  vortex  sheets,  which  offer 
further  evidence  for  this  non-uniqueness  of  the  Inviscid  flow.  With  these  symmetric 
separation  lines,  substantial  flow  asymmetry  is  still  possible  and  side  forces  similar 
in  magnitude  to  the  normal  force  are  predicted.  This  means  that  there  is  an  inviscid 
mechanism  which  is  capable  of  producing  the  observed  levels  of  side  force  with  no  appeal 
to  viscous  effects. 

Of  course,  the  real  separation  lines  are  not  symmetrically  placed,  so  the  third 
step  has  been  to  take  the  observed  separation  positions  from  an  experiment  and  supply 
them  to  the  inviscid  model  in  order  to  calculate  the  side  force.  As  indicated  above,  if 
the  solutions  are  of  the  first  family,  in  which  the  solution  reverts  to  its  symmetrical 
form  when  the  separation  lines  are  symmetric,  the  side  force  calculated  is  small,  much 
smaller  than  that  observed  in  the  experiment.  On  the  other  hand,  if  the  solutions  are 
taken  from  the  second  family,  which  remain  asymmetric  even  with  symmetric  separation 
lines,  the  calculated  side  force  agrees  quite  closely  with  the  measured  value. 

2  MATHEMATICAL  MODELS  OF  THE  FLOW 

The  simpler  of  the  models  used  is  an  extension  of  the  line-vortex  model  of  Bryson5. 
The  more  complex  model  is  a  development  of  the  inviscid  component  of  the  method  described 
in  Ref  3.  The  more  complex  model  used  is  described  first,  followed  by  the  simplifications 
that  reduce  it  to  an  extension  of  Bryson's  model. 

The  free  shear  layers  are  represented  by  vortex  sheets,  with  the  inner  turns  of  each 
vortex  sheet  replaced  by  an  isolated  line-vortex.  The  isolated  line-vortex  represents  the 
core,  and  is  joined  to  the  free  end  of  the  vortex  sheet  by  a  cut  as  indicated  for  the  star¬ 
board  vortex  in  Fig  1.  Boundary  conditions  are  applied  on  the  sheet,  such  that  it  is  a 
stream-surface  of  the  three-dimensional  flow  sustaining  no  pressure  jump.  The  conditions 
applied  to  each  line-vortex  and  cut  are  that  the  combination  is  force-free.  A  condition 
analogous  to  a  Kutta  condition  is  applied  at  the  separation  line.  The  vortex  sheet  model 
used  is  essentially  that  described  by  Smith6 ,  apart  from  differences  in  discretisation 
and  details  of  the  Kutta  condition.  As  in  Ref  3,  the  assumptions  of  slender-body  theory 
and  conical  flow  are  still  used,  but  the  flow  is,  of  course,  no  longer  symmetrical.  The 
slenderness  assumption  reduces  the  problem  to  the  solution  of  the  two-dimensional  Laplace 
equation  in  each  cross-flow  plane,  with  appropriate  boundary  conditions  on  the  vortex 
sheets  and  cone  surface.  The  assumption  that  the  flow  is  conical  means  that  we  need  only 
find  the  solution  in  one  cross-flow  plane.  Conformal  transformation  techniques  are  used 
to  simplify  the  solution  of  the  resulting  boundary-value  problem.  The  mappings  used  in 
Ref  3  relied  on  the  lateral  symmetry  of  the  flow  so  this  means  that  a  different  approach 
must  be  adopted  for  the  construction  of  the  complex  potential  describing  the  asymmetric 
flow  in  a  cross-flow  plane.  A  new  transformation  sequence  results. 
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The  velocity  field  ie  linear  in  its  components.  Fig  2  shows  the  decomposition  of 
this  field  into  the  contributions  arising  from: 

(a)  the  displacement  effect  of  the  growth  of  the  cone  cross-section  along  its  length 

(b)  the  angle  of  incidence  (or  incidence  plus  yaw)  of  the  cone 

(c)  the  isolated  vortices  representing  the  inner  parts  of  the  rolled  up  vortex 
sheet 

(d)  the  outer  part  of  the  vortex  sheets. 

Since  the  cross-section  of  the  cone  is  circular,  the  contributions  of  (a) ,  (b)  and  (c) 
to  the  velocity  field  can  be  written  down  at  once  as: 


If  we  describe  the  flow  by  this  velocity  field  and  apply  the  same  conditions  as  Brysons, 
ie  that  the  velocity  vector  at  the  cone  surface  at  each  separation  line-  lies  along  the 
separation  line,  and  that  the  total  force  on  each  combination  of  i  line-vortex  and  the 
cut  joining  it  to  its  separation  line  is  zero,  then  we  have  an  extension  of  Bryson's 
model  to  asymmetric  flow.  Details  may  be  found  in  Ref  4.  Some  results  from  using  this 
model  will  be  discussed  in  the  next  section. 

To  obtain  the  complete  model,  the  contribution  from  the  vortex  sheets  must  be  added 
to  the  expression  above.  For  reasons  discussed  in  Ref  3,  special  case  is  required  in 
representing  the  vortex  sheets.  Specifically,  the  curvature  of  the  vortex  sheet  at  its 
base  is,  in  general,  infinite  and  must  be  regularised  by  means  of  a  transformation  that 
maps  the  cross-section  of  the  cone  at  the  sheet  base  into  a  right-angled  concave  corner. 
Because  of  the  linear  nature  of  the  complex  potential,  each  vortex  sheet  can  be  considered 
separately  and  so  the  problem  reduces  to  the  canonical  one  of  finding  the  complex 
potential  due  to  a  single  vortex  sheet,  and  its  image,  leaving  a  circle  tangentially  and 
with  infinite  curvature.  This  is  solved  by  first  mapping  the  circular  contour  represent¬ 
ing  the  cross-section  of  the  cone  onto  the  imaginary  axis,  with  the  region  outside  this 
contour  being  mapped  onto  the  left  half-plane.  The  origin  is  then  shifted  to  the  base 
of  the  vortex  sheet  (the  separation  position)  and  the  axis  system  rotated.  A  folding 
(square  root)  transformation  is  then  applied  which  maps  the  imaginary  axis  onto  the 
positive  parts  of  the  real  and  imaginary  axes  and  the  right  half-plane  onto  the  first 
quadrant.  See  Fig  3.  The  construction  of  the  velocity  field  then  proceeds  as  explained  in 
Ref  3,  is  by  representing  the  sheet  as  a  collection  of  contiguous  circular  arcs,  with  con¬ 
tinuous  tangents  at  their  junctions,  in  the  final  transformed  plane,  Z5  .  Over  each  arc, 

a  distribution  of  vorticlty  is  specified  (which  is  also  continuous  across  the  junctions) 
such  that  the  velocity  field  for  the  sheet  may  be  written  in  closed  form.  In  effect,  we 
have  a  panel  method  to  represent  the  sheets. 

The  conditions  that  the  vortex  sheet  forms  a  three-dimensional  stream  surface  and 
that  the  pressure  is  continuous  across  it  are  applied  at  the  points  where  the  panels  join, 
and  at  the  mid-point  of  the  panel  nearest  the  core.  The  position  and  strength  of  a 
vortex  panel  are  specified  by  one  parameter  each,  so  that  these  conditions  are  just 
•officiant  to  determine  these  unknowns.  The  position  of  each  line-vortex  is  specified  by 
two  coordinates,  so  that  requiring  the  two  components  of  the  total  force  on  the  vortex  and 
the  cut  joining  it  to  the  free  end  of  the  sheet  be  zero  provides  equations  determining  the 
position  of  each  vortex.  The  final  unknowns  are  the  circulation  of  the  line-vortices. 

These  are  determined  by  the  Xutta  condition  that  the  velocity  vector  at  the  base  of  each 
vortex  sheet,  between  the  sheet  and  the  cone  surface,  lies  along  the  separation  line.  The 
resulting  sat  of  non-linear  algebraic  equations  is  solved  in  an  iterative  manner,  using  a 
multi-dimensional  analogue  of  the  Newton-Raphson  iteration  method.  In  this  method,  a 
goess  is  made  for  each  sheet  shape,  and  its  vorticity  distribution,  for  given  values  of 
the  incidence  parameter  and  separation  positions.  This  initial  approximation  is 
successively  Improved  by  the  Newton  method,  until  the  boundary  conditions  on  the  vortex 
sheets  are  satisfied  to  within  a  specified  error  tolerance.  Details  may  be  found  in 
Ref  3. 

3  RESULTS 

An  obvious  choice  for  an  initial  approximation  to  an  asymmetric  solution  is  one  of 
the  synaetric  solutions  found  previously3,  and  indeed,  the  first  solutions  for  asymmetric 
vortex  sheets  were  found  by  perturbing  the  symmetrical  solutions.  In  particular,  a 
gradually  increasing  asymmetry  in  the  separation  lines  was  introduced,  and  the  resulting 
solutions  found.  A  number  of  solutions  were  obtained,  covering  a  wide  range  of  separation 
positions.  No  attempt  was  made  to  use  any  particular  values  corresponding  to  experiment¬ 
ally  observed  separation  positions  at  this  stage  -  a  range  of  values  was  taken  in  the  hope 
of  covering  the  likely  range  of  experimental  results.  Of  particular  interest  is  the 
calculated  side-force  coefficient  Cy  with  asysmMtric  separation.  In  rig  4,  the  values 

of  the  ratio  Cy  to  the  normal  force  coefficient,  ,  are  shown  for  a  particular  value, 
a/4  -  3.5,  of  the  ratio  of  the  angle  of  incidence,  a  ,  to  the  semi-angle  of  the  cone,  4 . 
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The  abscissa  is  the  difference  in  angular  position  between  the  separation  lines  on  the 
port  and  starboard  sides,  so  that  the  origin  corresponds  to  zero  CY  with  symmetrical 

separation  lines.  The  two  curves  are  for  different  fixed  locations  of  the  port  separation 
lines.  The  most  significant  feature  is  that  the  value  of  C^/C^  is  small,  and  never 
exceeds  0.09. 

The  values  of  Cy/CN  measured  on  slender  pointed  bodies  in  general,  and  on  circular 

cones  in  particular,  depend  on  the  orientation  of  the  body  in  roll,  and,  in  some  experi¬ 
ments,  the  Reynolds  number.  The  largest  values  of  Cy/^jj  attained  as  the  body  is  rolled 

about  its  axis  are,  of  course,  of  the  greatest  interest.  For  circular  cones,  these  maximum 
values  depend  mainly  on  the  incidence  parameter,  a/i  .  For  the  value  of  this  parameter 
corresponding  to  the  results  shown  in  Fig  4  < 3 . 5 )  ,  the  maximum  values  measured  by 
Keener,  et  al1  exceed  0.5.  Since  the  range  of  separation  positions  covered  in  Fig  4  seems 
likely  to  include  any  occurring  in  a  real  flow,  it  is  clear  that  these  asymmetric  solutions 
do  not  correspond  to  the  observed  phenomena . 

If  the  solution  to  the  inviscid  flow  problem  were  uniquely  determined  by  the  body 
shape,  incidence  and  separation  line  positions,  that  would  be  the  end  of  the  matter. 

However,  for  such  non-linear  problems,  there  is  always  the  possibility  of  non-unique 
solutions,  although  finding  them  may  be  difficult.  The  vortex  sheet  model  is  quite  complex, 
and  carrying  out  an  extensive  survey  of  its  solutions  would  be  a  very  time  consuming 
process.  So,  the  simpler  line-vortex  method  was  used  to  search  for  other  possible 
solutions.  As  reported  in  Ref  4,  the  hunt  for  solutions  alternative  to  those  of  the  first 
family  was  successful,  and  a  new  family  of  solutions  was  found.  In  particular,  asymmetric 
solutions  were  found  for  separation  lines  placed  symmetrically  about  the  incidence  plane. 
These  asymmetrical  solutions  with  symmetric  separation  lines  typify  a  distinct,  second 
family  of  solutions;  the  solutions  obtained  by  Bryson  represent  the  first  family.  This 
classification  into  families  can  be  extended  to  solutions  with  asymmetric  separation  lines, 
in  an  obvious  way.  For  a  particular  value  of  the  incidence  parameter,  a  sequence  of 
neighbouring  solutions  can  be  found  (by  the  Newton-Raphson  procedure)  in  which  the  posi¬ 
tions  of  separation  are  varied.  If  the  members  of  this  sequence  with  symmetrical  separa¬ 
tion  lines  have  symmetric  vortices,  then  the  whole  sequence  of  solutions  belongs  to  the 
first  family.  If  not,  and  asymmetrical  vortices  are  still  found  when  the  separation  lines 
become  symmetrical,  then  the  sequence  belongs  to  the  second  family  of  solutions.  It  turns 
out  that  solutions  of  the  second  family,  defined  in  this  way,  give  rise  to  values  of 
Cy/Ctf  approaching  unity. 

Some  typical  results  are  shown  in  Fig  5,  where  the  locus  of  the  vortex  positions  for 
varying  incidence  and  fixed  (symmetrical)  separation  positions.  Two  sets  of  results  are 
plotted.  Those  labelled  1  are  members  of  the  first  family  of  solutions  and  are  symmetrical 
for  all  values  of  the  incidence  parameter,  because  the  separation  position  is  symmetrical. 
Those  labelled  2  are  examples  of  the  second  family,  where  the  vortices  are  not  symnetric. 

Fig  6  shows  the  corresponding  results  for  Cy/Cjj  calculated  from  the  second  family  of 

line-vortex  solutions.  As  can  be  seen,  values  similar  to  those  observed  in  experiments 
(in  this  case,  maximum  values  recorded  during  experiments  by  Keener,  et  al7)  ,  are  obtained. 

Eventually,  for  sufficiently  large  incidences,  the  results  of  the  line-vortex  model 
become  physically  unrealistic,  with  one  vortex  moving  a  very  long  way  from  the  cone.  Also, 
the  Kutta  condition  used  here  merely  enforces  a  crossflow  stagnation  point  at  the  separa¬ 
tion  line  position,  but  cannot  ensure  that  the  flow  is  away  from  the  body  at  this  point, 
ie  a  real  separation  line  behaviour.  In  fact,  at  large  incidences,  the  flow  is  found  to 
attach  at  these  points  instead,  which  is  another  aspect  of  the  physical  unreality.  Even  at 
values  of  the  incidence  parameter  which  are  of  practical  interest,  the  line-vortex  model 
is  unsatisfactory  for  several  reasons.  Most  fundamentally,  the  possibility  of  determining 
the  separation  line  in  any  rational  way  and  without  empiricism,  depends  on  a  detailed 
representation  of  the  flow  near  the  separation  line,  as  shown  in  Ref  3.  Even  in  the 
present  context,  where  we  are  treating  the  flow  as  inviscid,  with  prescribed  separation 
lines,  the  predictions  of  the  line  vortex  model  are  only  qualitatively  useful.  The  dis¬ 
crepancies  between  the  line-vortex  and  the  vortex  sheet  predictions  for  sharp-edged 
slander  wings  are  familiar;  they  are  greater  for  smooth  bodies. 

There  is  therefore  a  strong  incentive  to  seek  solutions  for  the  vortex  sheet  model 
which  correspond  to  the  second  family  of  solutions  of  the  line-vortex  model.  There  was 
also  the  hope  that  the  line-vortex  solutions  would  prove  useful  in  finding  the  vortex 
sheet  solutions.  This  eventually  proved  to  be  the  case  but  not  in  any  simple  way. 

Details  of  the  rather  complex  process  used  will  be  given  in  a  later  publication.  For 
the  present  we  confine  ourselves  to  an  examination  of  some  typical  results. 

Fig  7  shows  a  symmetrical  vortex  sheet  solution,  for  an  incidence  parameter  of  3, 
and  adjacent  to  it  an  asymmetric  solution  for  the  same  incidence  parameter  and  the  same 
symmetric  separation  positions.  Hence  this  asymmetric  solution  is  a  member  of  the  second 
family  of  solutions,  while  the  symmetric  solution  is  a  member  of  the  first  family.  The 
different  extent  of  the  outer  part  of  the  vortex  sheet  in  the  two  solutions  is  not  thought 
to  be  significant.  The  vortex  which  lies  closer  to  the  surface  is  vary  similar  in  shape 
to  the  symmetrical  vortices.  The  farther  vortex  is  very  different  with  a  long  length  of 
almost  straight  sheet,  rolling  up  at  its  free  edge,  rather  like  the  trailing  vortex  sheet 
from  a  slender  delta  wing  at  very  large  angles  of  yaw  (see  Jones1).  Also  shown  as  small 
circles  are  the  positions  of  the  vortices  in  the  corresponding  solution  of  the  line-vortex 
model,  ie  for  the  same  incidence  parameter  and  separation  line  locations.  It  is  apparent 
that  these  vortices  do  not  in  any  way  represent  the  core  position  of  the  vortex  sheet 
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solution  as  had  baan  hopad.  Further  examination  of  tha  vortex  shaat  solution  reveals  that 
the  bulk  of  the  circulation  resides  on  tha  portion  of  tha  shaat  represented  explicitly  - 
not  in  tha  core  itself.  This  is  particularly  marked  for  the  farther  vortex.  Consequently, 
in  going  from  the  line-vortex  solution  to  the  vortex  sheet,  a  major  re-distribution  of 
circulation  takes  place,  as  well  as  a  large  shift  in  core  position. 

Peake,  Owen  and  Johnson  report  some  experimental  results  for  the  asymnetric  flow  over 
a  circular  cone  in  Ref  9.  Included  in  their  results  are  vapour  screen  photographs  of  a 
cross-section  of  tha  flow,  showing  tha  two  primary  vortax  cores.  It  is  interesting  to  make 
a  qualitative  comparison  between  one  of  these  vapour  screen  photographs  (Fig  8)  and  a 
calculated  vortex  sheet  for  a  similar  configuration  (Fig  7) .  Note  that  the  observed 
asymmetric  flow  has  one  vortex  core  much  higher  than  tha  other,  and  both  cores  lie  close 
to  the  incidence  plane,  as  predicted  by  the  vortax  sheet  model. 

In  order  to  make  a  comparison  with  experimental  side-force  data,  we  choose  recent 
measurements  made  by  Mundell  at  RAE2.  He  has  appreciated  the  importance  of  rigid  model 
mounting  and  of  making  a  thorough  investigation  of  roll  angle  dependence  from  the  work  of 
Hunt  and  Dexter10  on  ogive-cylinders.  In  addition,  there  are  two  feature?  of  his  experi¬ 
ment  which  make  his  results  particularly  suitable  for  comparison  with  the  predictions  of 
the  present  theoretical  model.  The  wind  tunnel  model  (Fig  9)  consists  of  a  conical  nose 
blended  smoothly  into  a  cylindrical  afterbody.  The  model  is  constructed  in  two  portions 
and  the  forces  on  the  front  portion  are  measured  separately.  This  portion  is  conical  and 
sufficiently  far  forward  of  the  blending  section  for  the  flow  over  it  to  be  very  nearly 
conical,  even  at  the  low  speeds  of  his  experiment.  The  second  feature  is  that,  at  flow 
speeds  of  up  to  50  m/s,  the  boundary  layers  on  the  conical  nose  are  laminar  at  separation, 
so  that  the  separation  line  is  almost  straight.  Mundell  established  that  the  side  force 
for  a  given  Incidence  changed  sign  as  the  roll  angle  varied,  but  that  the  absolute  value 
of  the  side  force  was  almost  the  same  over  moat  of  the  roll  angle  range.  In  order  to 
reproduce  this  level  without  the  need  for  roll  angle  variation,  he  fixed  the  behaviour  of 
the  siu-  force  by  introducing  a  small  protuberance  at  the  apex  of  the  cone.  This  had  a 
negligible  effect  on  the  magnitude  of  the  side  force  and  also  made  it  possible  to  measure 
the  separation  line  positions  which  correspond  to  the  side-force  values,  using  an  oil  flow 
technique.  It  is  thus  possible  to  compare  Mundell* s  experimental  results  and  the  predic¬ 
tions  of  tha  vortax  sheet  model  in  some  detail.  He  shall  concentrate  on  a  comparison  of 
the  maxlsnim  side  force  at  a  particular  Incidence,  36  degrees.  For  this  case,  Mundell 
measured  separation  positions  at  0  degree  and  166  degrees.  It  is  Interesting  to  note  that 
this  asymmetry  in  separation  position  is  only  14  degrees.  These  separation  positions  were 
used  to  calculate  vortex-sheet  solutions  for  an  incidence  parameter  of  3.6  (Mundell *s  cone 
had  a  10  degree  semi-apex  angle) .  Two  solutions  were  found,  one  from  the  first  family, 
and  one  from  tha  second.  They  are  shown  in  Fig  10.  It  can  be  seen  that  the  solution  from 
tha  first  family  is  only  slightly  asysaatr.tc,  in  contrast  to  the  marked  asymmetry  of  the 
second  family  solution.  Comparison  with  the  solutions  of  Fig  7  shows  that  the  large  values 
of  the  incidence  parameter  and  the  earlier  separation  have  resulted  in  larger  vortices, 
as  would  be  expected.  At  the  same  time,  there  is  a  clear  'family*  resemblance  between  the 
solutions  for  the  two  cases,  despite  the  asysnetry  in  the  separation  position  In  Fig  10. 

A  furthar  comparison  of  tha  two  solution*  given  below,  in  tabular  form: 


First 

family 

Second 

family 

Expt. 

left 

BUI 

-13.86 

1  total 

right 

By;., 

14.02 

rcore^rtotal 

left 

0.80 

right 

0.60 

0.25 

v*2 

-1.53 

-19.73 

-17.6 

v*2 

36.81 

36.75 

27.0 

We  note  first  that  the  strengths  of  all  the  vortices  are  similar,  and,  in  particular,  that 
the  very  different  looking  vortices  of  the  second  family  have  almost  the  same  strength. 

The  next  line  of  the  table  gives  a  simple  comparison  of  the  way  in  which  the  circulation 

is  distributed  In  the  vortex,  r _ is  defined  as  the  circulation  of  the  isolated  vortex 

cor* 

plus  the  circulation  of  the  portion  of  the  sheet  which  lies  beyond  its  highest  point. 

There  Is  little  difference  in  tha  distribution  for  the  vortices  of  the  first  family 
solatia.::  but  in  the  second  family  we  see  that  tha  nearer  vortex,  on  the  left,  has  the 
bulk  of  its  circulation  in  the  core,  while  the  further  vortex,  on  the  right,  has  the  bulk 
of  the  circulation  in  the  sheet. 

The  force  coefficients  are  based  on  the  free-stream  kinetic  pressure  and  the  plan 
area  of  the  model  beck  to  its  split  line,  and  the  tabulated  quantities  are  further  divided 
by  the  square  of  the  semi -eagle  of  the  cone  (in  radians),  to  give  similarity  variables. 

As  expected,  the  side  force  from  the  first  family  Is  smaller  than  the  measured  value  by 
an  order  of  magnitude.  The  solution  from  the  second  family  gives  an  approximation  to  the 
measured  value  whioh  is  surprisingly  close,  considering  the  assumption  of  asiall  incidence 
which  is  inherent  in  slender-body  theory.  There  can  be  little  doubt  that  the  second 
family  corresponds  to  the  observed  flow.  The  solution  of  both  families  predict  much  the 
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same  level  of  normal  force,  which  Is  much  greater  than  that  measured.  This  is  a  reflec¬ 
tion  of  the  limitations  of  the  small  incidence  assumption  implicit  in  slender-body  theory. 
Figs  11  and  12  show  pressure  distributions  round  the  surface  of  the  cone  for  the  calculated 
solutions  shown  in  Fig  10.  The  Independent  variable,  9'  ,  is  the  angular  position 
measured  from  the  windward  generator.  The  ordinate  is  the  difference  between  the  local 
pressure  coefficient  and  its  value  at  the  windward  generator,  converted  to  a  similarity 
2 

form  by  division  by  6  .  In  the  solution  from  the  first  family,  Fig  11,  two  suction 

peaks  can  be  seen  on  each  side  of  the  cone.  The  first  is  related  to  the  maximum  of  the 

circumferential  velocity  from  windward  to  leeward  which  occurs  near  the  equatorial  position. 
The  pressure  rise  following  this  suction  peak  leads  to  the  flow  separation,  following  which 
the  pressure  is  initially  almost  uniform.  The  second,  much  larger,  suction  peak  follows, 
corresponding  to  a  maximum  in  the  circumferential  velocity  from  leeward  to  windward  under 
the  vortex  core.  In  a  real  flow  the  second  peak  would  be  modified  by  the  effects  of 

secondary  separation  in  a  similar  way  to  that  observed  on  delta  wings.  In  this  solution 

from  the  first  family,  the  dominant  contribution  to  the  side  force  comes  from  the  difference 
between  the  first  peaks  on  the  two  sides  of  the  body.  The  contribution  from  the  second 
peaks  is  smaller,  because  the  lateral  projection  of  the  surface  element  on  which  it  acts  is 
much  smaller. 

In  the  solution  from  the  second  family,  shown  in  Fig  12,  the  most  striking  feature  is 
the  absence  of  any  suction  peak  related  to  the  farther  vortex  core  shown  in  Fig  10.  At  the 
same  time,  the  suction  peak  of  the  nearer  core  is  much  enhanced,  as  is  the  suction  peak 
preceding  separation  on  the  side  where  the  core  is  nearer.  The  side  force  is  consequently 
much  larger  for  the  second  family,  in  agreement  with  the  overall  forces  presented  earlier. 
Another  feature  of  the  increased  asysmetry  of  the  solution  of  the  second  family  is  the 
displacement  of  the  leeward  attachment  line  by  some  20  degrees  from  the  leeward  generator. 
The  positions  of  the  other  attachment  and  separation  lines  are  also  indicated  in  Figs  11 
and  12. 

The  results  presented  here  are  very  much  those  of  'work  in  progress',  and  a  more 
extensive  survey  of  the  second  family  of  vortex  sheet  solutions  is  needed.  In  particular, 
it  should  be  possible  to  discover  for  what  range  of  incidence  parameter  and  separation 
positions  the  second  family  of  solutions  exists,  and  whether  they  can  be  related  to  the 
first  family  in  any  way. 

Such  a  survey  has  already  been  conducted  with  the  simpler  line-vortex  model.  Here, 
the  second  family  of  solutions  are  found  to  arise  from  the  first  family  by  a  process  of 
bifurcation.  This  is  illustrated  in  Fig  13,  where  the  coordinates  of  the  vortex  are 
shown  as  functions  of  the  incidence  parameter,  the  separation  positions  being  fixed 
syamietrically.  The  curves  labelled  'symmetric'  represent  the  solution  for  the  first 
family.  For  symmetric  separations,  these  are  Bryson's  solutions,  which  start,  with  a 
vertical  tangent,  where  a/a  «  1.5  cosec  8  *  1.81  .  At  a  value  of  the  incidence  parameter 
just  above  2,  at  the  point  labelled  B  ,  the  asymmetric  solutions  depart  from  the 
synnetrlc  solution  along  a  vertical  tangent,  in  opposite  directions,  when  attempting  to 
calculate  the  solutions,  it  is  important  to  consider  that  the  determinant  of  the  Jacobian 
matrix,  on  which  the  Newton  scheme  depends,  vanishes  at  B  .  We  see  from  Fig  13  that  the 
second  family  of  solutions  for  the  line-vortex  model  does  not  exist  for  all  combinations 
of  separation  position  and  incidence  parameter.  Fig  14  shows,  for  symmetrically  placed 
separation  lines,  (81  +■  ,  -  180  degrees)  ,  two  boundaries  in  the  plane  of  incidence 

parameter  a/6  ,  and  separation  position,  9 ^  .  The  first  is  the  boundary,  found  by 

Bryson5  below  which  there  are  no  solutions  for  the  first  family  with  the  vortices  near 
the  surface  of  the  cone.  The  second  is  the  bifurcation  boundary,  below  which  no 
asymmetric,  second  family  solutions  are  found.  It  is  plausible  that  similar  boundaries 
for  the  existence  of  solutions  of  the  vortex  sheet  model  falling  into  the  first  and 
second  families  exist,  but  more  work  is  needed  to  establish  them.  We  note  that  symmetric 
vortex  sheet  solutions  are  found  in  Ref  3  at  values  of  the  incidence  parameter  which  are 
below  the  lower  limit  for  solutions  of  the  line-vortex  model. 

4  CONCLUSIONS 

It  is  important  to  recall  the  limitations  and  approximations  of  the  present  approach 
to  the  problem  of  out-of-plane  forces  on  bodies.  In  the  first  place,  attention  has  been 
concentrated  on  the  maximum  levels  of  side  force  which  occur.  Secondly,  consideration  has 
been  confined  to  circumstances  in  which  separation  along  the  body  sides  extends  all  the 
way  forward  to  the  pointed  apex  of  the  body,  in  other  words,  for  angles  of  incidence  which 
are  large  in  terms  of  the  semi-angle  of  the  apex.  Out-of-plane  forces  also  arise  in  other 
circisnstances,  for  instance  through  separation  on  cylindrical  afterbodies.  For  flows 
which  do  separate  from  the  apex,  the  present  treatment  is  confined  to  circular  cones, 
though  it  could  readily  be  extended  both  to  cones  of  other  cross-sections  and  to  non- 
conlcal  bodies.  However,  the  limitations  inherent  in  slender-body  theory  would  have  to  be 
borne  in  mind  in  any  attempt  to  make  use  of  its  quantitative  predictions  in  aerodynamic 
design. 

What  is  claimed  for  the  present  treatment  is  that  it  reveals  an  inviscid  mechanism  - 
naswly  the  existence  of  multiple  inviscid  solutions  to  a  simple,  well-defined  problem  - 
which  can  give  rise  to  flows  and  forces  which  ars  close  to  those  observed,  the  results 
suggest  that  this  mechanism  la  responsible  for  tha  occurrence  of  out-of-plane  force  on  the 
whole  class  of  pointed  bodies  under  conditions  in  which  separation  extends  to  the  apex, 
though  tha  downstream  developsMnt  of  the  flow  will  inevitably  be  much  more  complicated 
than  the  conical  structure  assumed  here. 
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